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Abstract 

We study the lepton flavor violation (LFV) in tau decays in the framework of the super symmetric 
seesaw mechanism with nonholomorphic terms for the lepton sector at a large tan (3. In particular, 
we analyze two new decay modes r — > £/o(980) and r — ► £K + K~ arising from the scalar boson 
exchanges contrast to r — » trf^ from the pseudoscalar ones. We find that the decay branching 
ratios of the two new modes could be not only as large as the current upper limits of O(10 -7 ), but 
also larger than those of r — > irff'. Experimental searches for the two modes are important for the 
LFV induced by the scalar-mediated mechanism. In addition, we show that the decay branching 
ratios of r — > lu + u~ are related to those of r — > In and r — > ^/o(980). 
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In the standard model (SM), since the neutrinos are regarded as massless particles, the 
processes associated with lepton flavors are always conserved. Inspired by the discoveries of 



nonzero neutrino masses 111 I2JI , it has been studied enormously how to generate the neutrino 
masses which are less than a few eV. By supplementing with singlet right-handed Majorana 
neutrinos with masses Mr required to be around the scale of unified theory, it is found 
that the seesaw mechanism is one of natural ways [3] to obtain the small neutrino masses. 
Accordingly, in non-SUSY models, it is easy to understand that the effects of the lepton 
flavor violation (LFV) are suppressed by l/M R . However, in models with SUSY, due to the 
nondiagonal neutrino mass matrix, the flavor conservation in the slepton sector at the unified 
scale will be violated at the Mr scale via renormalization P,^, f|. The flavor violating effects 
could propagate to the electroweak scale so that instead of l/M R , the suppression of the 
LFV could be 1/Msusy with Msusy ~ 0(TeV) being the typical mass of the SUSY particle. 
Consequently, the lepton flavor violating processes, such as £~ — > £j'j and £~[ — > £ji££k, 
become detectible at the low energy scale. The LFV has been extensively studied in the 
literature. For examples , r — > \ir\, B — ► (e, fi)r and the \i — (e, r) conversions can be found 



in Refs. 



Ref. (l3|. 



aaayy 



121 ]. while that to the detection of the LFV in colliders is given in 



In the large tan /3 r egio n, it has been pointed out that the nonholomorphic Yukawa 



interactions 1^ 



15. 



p r egio 

HQ 



play very important roles for flavor changing neutral currents 



m 



(FCNCs) in the quark sector. In the SUSY-seesaw model, the nonholomorphic terms 
the lepton sector naturally induce the LFV due to the Higgs couplings. It has been shown 
that the contribution to the decay of r — > 3/x from the Higgs- mediated LFV at large tan (3 
could be much larger than that from r — > /ry — > 6|, [l8|. Recently, the experimental 

limits on the radiative decays of r — ► £7 (I = e, /1) have been improved from O(10 -6 ) ^| 



to O(10 7 ) 20j, |2l|. Moreover, the sensitivity of probing the LFV in r decays with single 



pseudoscalar (P) or vector (V) and double mesons in the final states, i.e., t — > £(P, V) and 



t — > £PP, have also reached O(10 
t — > 



221 ] . In this paper, we will simultaneously analyze 
and r — > £X, where X are rj^'\ (J), / (980), cr(600) and K + K~ , respectively, in 



the Higgs-mediated mechanism. In particular, we would like to check whether it is possible 
to have large rates for the processes beside the mode of r — > 3fi. Note that the decays 
of r — » £S with S = / (980) and a(600) and r — > £((fi, K + K~) have not been explored 
previously based on the Higgs-mediated mechanism in the literature, while r — > £P have 
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been studied in Refs. 

We start with the Higgs-mediated mechanism. It is known that by the induced slepton 
flavor mixing, the effective Lagrangianjvith induced nonholomophic terms for the Higgs 
bosons coupling to leptons is given by 



£ eff = E R ^ [5 l3 H° d + (e 1 5 ij + e^) H° u *] E Lj + h.c. , 
= E R M°E L + h.c. , 



(1) 



where Y denotes the diagonalized Yukawa matrix of leptons, 7y = (Am|)y/mg_and ei(2) is 
related to the induced lepton flavor conserving (violating) effect, expressed by 



e a ~ —fxM, 2/ 1 (M 1 2 , mf T , mfj - f x {M(, /T, mf T ) + 2f 1 (M(, mj r , mjj 
an 

a 2..u ,2 r ; ,.2 m 2 ,,,2 1/-2, o-ff.,2 ,,,2 ,,,2 i f2 



+ 



/i(// 2 ,mf r +2/ 1 (/i 2 ,ml f ,M|) 

2 2 2 \ _ 

l e L i m fL' m fRi J*\r i ,,v e. L 
>M 2 m 2 / 2 ( M 2 , mi , m 2 L , M 2 ) + 2/ 2 (/i 2 , m 2 f , m^, M 2 ) 



(2) 



where M x 2 are the masses of gauginos from the soft SUSY breaking terms, \x stands for the 
mixing of H u and Hd, 

xy\n(x/y) + yzln(y/z) + zxln(z/x) 



fi(x,y,z) 
f 2 (w,x,y,z) 



(x - y){y - z){z - x) 

MilllMJ 

— cyclic, 



(w — x)(w — y)(w — z) 

and a 1(2) = 9xi 2 \/^ with <?i( 2 ) corresponding to the gauge coupling of the £/(l) (577(2)) 
symmetry. Due to the nonholomorphic term e 2 Eij, the lepton mass matrix is not diagonal 
anymore. Consequently, after rediagonalizing the lepton mass matrix, the lepton flavor 
changing neutral interactions through the Higgs bosons appear. Since the nonholomorphic 
terms are expected to be much less than unity, to obtain the LFV, we take the unitary 
matrices used for diagonalizing lepton mass matrix to be Ul(r) ~ 1 + ^l(r) as a leading 
expansion of e 2 i%, where Al(r) are 3x3 matrices. From Eq. ((THj) and (Am|)y = (Am|)jj, 
we may set Al = A R = A. Hence, the diagonal mass matrix in Eq. (P) could be obtained 
by 



UM$lP « (1 + A)M°(1 - A) = M; 



din 
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where Mf ia is the physical mass matrix of the lepton with the diagonal elements being 
(Mf ia )u = (m e , rrifj,, m T ). At the leading order, we get 



In terms of the physical mass eigenstates of the Higgs bosons, represented by (24 1 



ReH® = v d + -j= [cos aH° — sin ah ] , ReH® = v u + — ^= [sin ai/ + cos , 
ImH% = [cos(3G° - sin (3 A ] , Jm#° = -j= [sin/3G° + cos/M°] , 

where a is the mixing angle of the two CP-even neutral scalars, the interactions for the LFV 
via the Higgs-mediated mechanism are expressed by 

H% j = (v^GV) 1/2 ^4f [sin(a - /3)H° + cos(a - (3)h° - iA°] + h.c. (3) 

COS l) 

with ma is the mass of the ith flavor lepton and Cy = €%lijj (1 + (ei + tan/3) 2 . 

From Eq. we see that the decays of r — > £P only pick up the contributions from the 
pseudoscalar boson A , while t ^ £S and r — > £PP are governed by both scalar bosons 
P° and h° due to the parity properties. In our following analysis, we only concentrate on 
the processes associated with the productions of ss and pairs to avoid small Higgs 

couplings. We choose the decays of r — ► £X with X = yU + /i~, 77 ^, /o(980)(cr(600)) and 
K + K~ as the representative modes. For r — * £fi + fi~, the formalisms for the decay rates 
dictated by scalar and pseudoscalar bosons are given by 



G 2 m 2 ml\C Ti \ 2 



cs sc \ 2 / sin /3 N 
m 2 m 2 H ) V m A 



(4) 



3 • 2% 3 cos 6 (3 

where q = 3/2 and 1 with £ — /i and e, cs = cos(a — 0) sin a and sc = sin(o; — (3) cos a, 
respectively. To study the production of 77 we adopt the quark- flavor scheme, defined by 



cos </> — sin 
sin 6 cos 




(5) 



where r\ q = (uu + dd)/^/2 and i] s = ss. From (0| ^'7^759' = f Vq ,P^, the mass of r) q ^ 
can be expressed by m 2 qq = ^(0\m u u^5U+m c [d'y 5 d\ri g ) {m 2 s = j-{Q\m s s^s\r] s )). If we neglect 
the r\ q contribution due to small m u ^, the decay rates for r — > £7] can be written as 

F(r - ft,) - tan« „ (sin tt < ) 7l - ^) ' . (6) 

647T \ J \ m l T ) 



Similarly, the rate for r — > irj' is given by 



r(r — > £77) \ 1 — rn%jm\ 

For r — >■ £/o(980)(cr(600)) decays, although the quark contents of /o(980) and <r(600) are still 
uncertain, we adopt two quark contents to describe the states. In terms of the notations in 



Refs. 
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271 ]. the isoscalar states / (980) and cr(600) are described by |/ (980)) = cos6*|ss) + 
sin#|rm) and |c(600)) = — sin#|ss) + cos9\nn) with nn = [uu + dd)/y/2 and 9 being the 
mixing angle. The decay constants are defined as 

(/ *M0> = m f J s fo , (a s \ss\0) = mj s a , (8) 

where Jq and a s represent the ss component in / (980) and cr(600), respectively. As a result, 
the decay rates of r — > £/ (980) are given by 

^ - « - G i^k ( (5 - 1) 2 - S) 2 • (9) 

On the other hand, the rates for r — > £a(Q00) can be obtained by 

r(r -> £a(600)) / mJgtang V / 1 - m£/m? V 

r(r - £/ (980)) " ^ royb /. J 1^1 -m\l ml) ' 1 J 

For the three-body decays of r — > £K + K~ , the associated hadronic effects are much more 
complicated and unclear. Nevertheless, the uncertainties could be fixed by the S decays, 
such as B — > KKK. The related form factor including resonant and nonresonant effects is 
defined by j3] 

£s S^KK 

<a-+( Pi )a-( P2 )mo) , tr-vh = e m r^ Mnsrs + f r, (id 

where S 1 stands for the possible scalar meson state, ms/J = (*S'| ss |0) , g s ~^ KK denotes the 
strong coupling for S — > KiT, and 

12\ -1 



/f R = i(3^, + 2^)+^(ln|!) 



1(2) 




q 2 g 4 

with u = (m| - m 2 )/(m s - m d ), a;J = -3.26 GeV 2 , x\ = 5.02 GeV 2 , x\ = 0.47 GeV 2 
and x\ = 0. It is found that only / (980) and / (1530) have the largest couplings to the 



28| . the factorization approach in Ref. 



KK pair 291. Note that in calculating B -> KKK 

n 

[30( has been used. In our numerical estimations, we will only consider these two scalar 
contributions. The differential decay rates as a function of the invariant mass in the KK 
system are given by 



dT(r -> £K + K~ 
dQ 2 



x 



G 2 F m\\C Tl \ 2 
2% 3 cos 6 (3 



m s f t 



K+K- 



CS 



mi 



sc 



m 



li 



Am 2 K 



1/2 



(13) 



m* j \ Q 2 

From Eqs. (jlj), © and 0, it is interesting to see that the various decay rates mediated by 
the Higgs bosons have the relationship 

Tfr -+lry) | Tjj -+ £/ (980))" 



r(r -f 



Qm 2 m^ 

3 • 2 5 7T 2 



a 



c 



/o 



(14) 



2\2 



where = (sin 2 j3 sin <fif s m 2 s /2) 2 (l — m 2 /m 2 )) 2 and C/ = (m s mf fj o cos #) 2 (1 — m 2 o /m; w 

We now consider the radiative modes of r — > £7. At the large tan (3 scenario, the dominant 
contributions to the decays are illustrated in Fig. ^ To simplify the estimations, we use the 
mass insertion method to formulate the decay amplitudes. The induced LFVs in the slepton 
mass matrix can be approximately written as 



amp ltud 

fly, Ei 



mi, 



1 



(6m 2 YX + ^A)..\n 



Mr, 



(15) 



■v (4tt) 2 v- ' Hj- y Mliy 

where m , Y u and A v denote the typical initial soft SUSY-breaking mass of the slepton, the 
neutrino Yukawa couplings and the trilinear soft SUSY-breaking effects, respectively, at the 
unified scale of Mjj. From Fig. and Eq. (|15|). the effective interactions for r — > £7 are 



, - x- - .. 



, - "*~ - - s 



X . X 



H" H 



TI-' 1 



X . X 



(a) 



HI HZ KG 
(b) 



FIG. 1: The Feynman diagrams for r — ► ^7 with a large tan/3. The corsses represent the various 
mixing effects. 



given by 



T 



G p 

71 



em T e^{k)t{p - k)ia^A R {\ + 75 )r(p) 



(16) 
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where 



.4 



R 



M 2 jJ. m 



w 



(4vr) 2 M| 



tan/3(Am|) T< ^ G £ 



S=£,u e 



(17) 



Gi 



1 — tan 2 9w 



2 



m 



2 

m~ 



TL 



2 

m~ 



m 



fc{%i>e) fc(%i> T ) 



2 



/n(z) 



(l — a; 2 + 2x In x) , f c (x) 



2(1 



(3-4x + x 2 + 21nx) , (18) 



with xs = M^/rrig |5|. Here, we have set the masses of higgisions and gauginos to be the 



same, denoted as M 2 . Subsequently, the decay rates of r — > £j are given by 



-GlmllA 



ft 



(19) 



The diagrams in Fig. ^ can also induce r — > £fi + fi , r — > £0 and r — > £K + K when the 
photon is off-shell. From Eq. ()16j) . it is easy to estimate the ratios of branching ratios (BRs) 
to be Q 



BR(r -> £X 



7/7 



10" 



(X, 



(20) 



BR(r^£ 7 ) \ 7T 

Note that it is impossible to produce modes with X 7 being a single pseudoscalar or scalar 
by the dipole operators in Eq. (fTT)J) . In our estimations for the modes with X 7 = <fi and 
K + K~ , we have used the hadronic matrix elements defined by (O|g7 M g|0) = im^fye^ik) and 



1 (77^(7 1 if + fa (p 2 )) = G°i — P2)Fq C+K (Q 2 ), with the form factors given in Refs. [32 



and 



(28j, respectively. It is clear that from the current limits on BR(t — > £7), BRij — > OT 7 ) 7 
are too small to be observed. We remark that other loop contributions to the decays, such 
as those from box diagrams, are expected to be small due to the light fermion final states. 

For the numerical estimations on r — > ^7 and r — > £X, we assume that M\ ~ M 2 ~ mo ~ 
/i ~ mi ~ mf to simplify our discussions. Consequently, Eqs. and (fTTj) become 



3a f 



47rsin 2 (2^) 



?2 



167T V 3 cos 2 6 



+ 



w 



sin 2 W 



-4, 



1 mi (Am^-Jrf 



tan/3(l + tan 2 9w), 



6(4-7r) 2 m? m^ 
respectively. If we regard A'A in Eq. (|15J) as (A^A)^ ~ m, 



(21) 



m,Q0(l), we get 

(Am 2 -) r£ /mg 8/(4vr) 2 ln(M u /M i? ). Thus, we find that C rl are insensitive to the SUSY 
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breaking scale and the decays of r — > £7 and r — > £X are only sensitive to the masses 
of the slepton and Higgs bosons, respectively. In calculating the numerical values, we set 
Gu(R) — 10 19 (10 14 ) GeV and tan (3 = 60. Other parameters in various modes are taken to be 



as follows: <\> = 39°, f s = 0.17 GeV and m ss = 0.69 GeV for r -> V } 25]; 9 = 30°, m s = 0.15 
GeV and /* ~ f s fo = 0.33 GeV (2T | for r -> £(/ (980), <r(600)); u = 2.87 GeV, k = -10.4 
/ /o(980) = 0.33 GeV, ghm^KK = L5Q GeV; g f (wao)-+KK = 3 18 GeV 



GeV 4 , // (i53o) 
|28|, r /o(980 ) = 



80 MeV and T 



/o(1530) 



1.16 GeV |29| for r -> tK + K~. For simplicity, we 



do not distinguish the difference between (Y^Y^)^ and (Y^Y^)^, i.e., (Am|) re = (Am|) T/J . 
In Fig. El we present the BRs for r — > £7 as a function of the slepton mass. In comparison 




7.2 1.5 

m~ (TeV) 



FIG. 2: Branching ratios (in units of 10 7 ) for r — > ^7 as a function of the stau mass. 



with the BELLE and BABAR results of Ei^r -> /ry) < 3.1 x 10~ 7 



21(, we see clearly that mf > 1 TeV is favorable. The BRs of r 



and 0.68 x 10~ 7 
£i] as a function 



of the pseudoscalar mass are displayed in Fig. 0(a). From Eq. (JJJ), we have BR(r — > 
&/) = 0.93BR(t -> £77). The BRs of r -> £/„(980) and r -> fiY+iY" as a function of 
M# = (cs/m\ — sc/m 2 H )~ 1 / 2 are shown in Figs. 0(b) and 0(c), respectively. In terms of 
Eq. (JIOD, we get Ei^r -> £a(600)) = 0.2BR(r -> £/„(980)). In addition, from Eq. (m, 
we obtain that £i2(r -> V>~) - °- 33 [ 5i? ( r ^ £ v) + l-GBR(r -> £/ (980)]. Clearly, all 
r — ► modes except r — > £er(600) are suitable to search for the LFV. Finally, it is worth 
mentioning that if we take the decoupling limit, i.e. m# « and a — > (3 — 7r/2 [24j, leading 
to = m ff , we get T(r -> £/ (980)) : T(r -»■ : T(r -»■ £77) « 1.3 : 0.36c* : 1. 

In summary, we have studied the lepton flavor violating r decays through the Higgs- 
mediated mechanism with the nonholomorphic terms from the couplings between the 
Higgs bosons and leptons at the large tan/3. By assuming that all masses associated with 
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FIG. 3: Branching ratios (in units of 1(T 7 ) for (a) r £r] and (b)[(c)] r -► ^/ (980)[if as 
functions of the pseudoscalar and scalar Higgs masses, respectively. 

SUSY-breaking are the same, we have demonstrated that BRs of r — > £7 only depend on 
the stau mass. In the Higgs-mediated mechanism, we have shown that the BRs of the new 
proposed decays of r — * ij " (980) and r — > IK + K~ arising from the scalar exchanges can 
be as large as the upper limits 0(1CT 7 ) of the current data and, moreover, they can be 
larger than those of r — > £rj from pseudoscalar exchanges. We have also pointed out that 
t — > £fi + fi~ are related with t —> £r] and r — > ^/ (980). It is clear that future experimental 
searches for the LFV in the leptonic and semileptonic tau flavor violating decays are 
important for us to identify the Higgs-mediated mechanism. 
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